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We study the charge transfer between a quasi-free-standing monolayer graphene, produced by 
hydrogen intercalation, and surface acceptor states. We consider two models of acceptor density of 
states to explain the high hole densities observed in graphene and find the density responsivity to 
the gate voltage. By studying magneto-oscillations of the carrier density we provide an experimental 
way to determine the relevant model. 
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Among numerious ways of graphene fabrication [TH3] , 
one of the promising methods for the top-down manu- 
facturing of electronic devices consists in the graphitiza- 
tion of Si-terminated surface of silicon carbide. Graphene 
produced by such graphitization resides on an insulat- 
ing (6-\/3 X 6-\/3) R30° carbon buffer layer, which is eas- 
ily spoiled by vacancies and, typically, leads to high 
graphene doping :[3]. Without special growth proto- 
cols, monolayer graphene on Si face of SiC appears 
to be highly n-doped, with electron density at about 
rie = 1 • 10^^ cm~^ level [51 [S]- Also, charged sur- 
face donors induce Coulomb scattering, which limits the 
mobility of electrons in such a material. To improve 
the transport qualities of graphene and decouple it from 
the substrate, it has been proposed to use hydrogen in- 
tercalation [7]. Hydrogen breaks Si-C bonds and con- 
verts the buffer layer into a quasi-free-standing mono- 
layer graphene (QFMLG), which is usually p-doped [8] 
due to electron transfer from graphene to acceptor states 
in the H-terminated surface of SiC. 



In this work we study the charge transfer between 
graphene and SiC in hydrogen-intercalated epitaxial 
graphene, including the magneto-oscillations of the hole 
density in QFMLG. We consider two limiting models for 
the charge transfer, which are different by the form of the 
density of states of surface acceptors, 7(e). In the model 
I we assume that all acceptor states are due to occa- 
sional unsaturated Si bonds: vacancies in the hydrogen 
layer with the density n„, energy Ei and narrow spec- 
tral density 7/(e) = nt,exp(-(e - EiY / /!^^)/{^jTi/\) ^■ 
nyS{e—Ej), A ^ Eg — Ej, where Eq is the work function 
of graphene (Fig. Ill inset I) . In the model II we assume 
a uniform density of states of acceptor levels 7(e) = 70 
filled up to the energy En (Fig. Ill inset II). For both 
models wc find the hole density in graphene and the re- 
sponsivity factor describing the effectiveness of QFMLG 
carrier density control by a gate voltage in a field-effect 
transistor. We also analyse the oscillations in the carrier 
density dependence on a magnetic field, and show that 
this can be used to distinguish experimentally between 
the two limiting doping models. 



The electron transfer from QFMLG to surface acceptor 
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FIG. 1: a) The dependence of the hole density n in QFMLG 
on the acceptor density (see n„ axis for the model I and 70 
axis for the model II). The insets show charge distribution 
between graphene and acceptor states for both models, b) 
Responsivity factor r dependence on riv The following pa- 
rameters were used for both plots: d = 0.3 nm, Aj — 0.6 eV, 
A = 5 meV, An = 0.7 eV. 



states is described by the following equation, 

e(n) 
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where n is the density of holes in graphene, Ug = CVg/e 
{Vg is gate voltage), d is the distance between SiC and 
QFMLG, the Fermi energy (relative to the graphene 
Dirac point) is ep^n) — —hvy^nn, and Emm = ^oo for 
the model I and Emin — En for the model II. Note 
that integral in Eq. T] is taken over the electron (rather 
than hole) energy. Each of the models is characterized 
by one energy parameter: A = Eg — Ej for model I and 
A = Eg- Eh for model II. 

The value of the carrier density for non-gated struc- 
tures can be obtained by solving Eq. (fTl) for Ug — 0. 
In Fig. [l^ we illustrate the hole density dependence in 
graphene on the amount of acceptors on hydrogenated 
SiC. For the model I, it is 



4^2 



n/ = mm [n.y, nf,J 



rib 



{Hv^/tt + Ji 



h^v'^i 



eo 



(2) 

At a small n„, all acceptor levels are occupied and nj = 
Uy. As acceptor density n^ increases, graphene doping 
saturates at n/ sa 1 • 10^'^ cm^^. Model II also shows 
a saturation of the carrier density, but with a smoother 
crossover, 



4^2 



nil = 



hv\ 



tfiv'^TT + AA 






(3) 



The effectiveness of using QFMLG transistors can be 
characterized by the responsivity factor, r — dn/dug. 
Then, r = 1 corresponds to the regime of the effective 
transistor operation, while r <?C 1 indicates that it is dif- 
ficult to change carrier density in graphene. The respon- 
sivity factors for the models I and II, 
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are compared graphically in Fig. [TId. 

Due to the charge transfer between graphene and ac- 
ceptors in the hydrogen layer on SiC, carrier density 
of QFMLG change upon the variation of the magnetic 
field B. For high acceptor densities [n^ ^ 10^^ cm"^^ 
7o ^ 10^^ cm~2eV~^), Fig. ^ shows the oscillations in 
n{B) dependence, which appear due to a non-zero density 
of acceptor states at the Fermi level. Similar oscillations 
were found in graphene on a non-hydrogenated SiC sur- 
face [5] , and used to explain the presence of a wide v = 2 
plateau in QHE. Both models give close results for the 
n{B) dependence, which reveals two different regimes, 
(a) The regions with a negative slope indicating the pin- 
ning of the Fermi level to one of the Landau levels in 
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FIG. 2: a) The magneto-oscillations of the hole density n in 
QFMLG in the regimes of large acceptor density, b) n{B) 
dependence at small acceptor density. The energy structure 
in different regimes is shown in the insets A-C (model II) and 
inset D (model I). The parameters used here were chosen for 
illistration purposes only. 



graphene, E]\j — —V2ehv'^BN. To evaluate n{B) depen- 
dence in this regime we solved Eq. (Il]), with ep = Em- 
(b) The regions with a positive slope corresponding to 
the pinning of filling factor v = nh/{eB) = AN + 2, for 
which the Fermi level in the system lies in the acceptor 
band, between Landau levels. 

The difference between the density oscillations in 
graphene expected on the basis of the models I and 
II is quite pronounced in structures with low acceptor 
densities {riy < 10^^ cin~^, 70 



< 



as 



shown on Fig. [2b. In this regime no oscillations are ob- 
served for model I, in contrast to the model II, where 
pinning of i^ = 2 takes place over magnetic field inter- 
val of several Tesla. In the latter case, one can iden- 
tify the following charge transfer regimes illustrated us- 
ing sketches in Fig. [2b. At high magnetic fields {B > 
Bi — hA/2e/{'jQ^ + e^d/eo)) the carrier density remains 
constant, since the 0th Landau level (LL) is partially oc- 
cupied and the Fermi level in graphene is pinned to the 
Dirac point (inset A in Fig. [2b) . At lower magnetic fields 



{2A^/{^fdn^+^Jehv^^ + 46^/(70"^ + e'^d/eQ)/hf <B< 
Bi) 0th LL becomes completely unoccupied and the 
Fermi level sticks to acceptor levels between 0th and 1st 
LLs (Fig. [2J3, inset B). This regime, observed within a 
window of several T, is responsible for widening the v — 2 
QHE plateau. The region C in Fig. [2]d refers to the pin- 
ning of the Fermi level to the 1st LL Ei . 

In conclusion, we have considered two models of the 
charge transfer between hydrogen intercalated graphene 
and acceptors in SiC surface. The first model describes 
the case when all acceptor states are created by unsat- 
urated Si bonds with energies from a narrow window. 
The opposite limit of a wide acceptor energy distribu- 
tion is covered by the second model. Both models pre- 
dict a saturation behavior of carrier density in graphene 



as the acceptor density increases, however, the magnetic 
field dependence of charge transfer in them is very dif- 
ferent, especially in structures with a low initial doping 
n < 10^^ cm~^ by holes. As shown in Fig. l2p, the model 
II features a wide (several T) i/ = 2 QHE plateau, while 
the model I reveals no ji{B) dependence. One of the ways 
to detect pinning of filling factor ;/ = 2 is to measure the 
activation energy T* and/or breakdown current of QHE 
i9j . The width of the peak in the dependence of T* and 
breakdown current on the magnetic field in the vicinity 
oi ly = 2 QHE plateau determines the region of filling fac- 
tor pinning. Alternatively, one can determine the filling 
factor by looking at the Landau level occupancy using 
scanning tunneling spectroscopy as in Ref. [10]. 
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